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We calculate the QCD equation of state for temperatures corresponding to the transition region 
with physical mass values for two degenerate light quark flavors and a strange quark using an im- 
proved staggered fermion action (p4-action) on lattices with temporal extent N T — 8. We compare 
our results with previous calculations performed at twice larger values of the light quark masses as 
well as with results obtained from a resonance gas model calculation. We also discuss the deconfining 
and chiral aspects of the QCD transition in terms of renormalized Polyakov loop, strangeness fluc- 
tuations and subtracted chiral condensate. We show that compared to the calculations performed 
at twice larger value of the light quark mass the transition region shifts by about 5MeV toward 
smaller temperatures. 
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I. INTRODUCTION 

First calculations of the equation of state (EoS) of hot 
strongly interacting matter date back to the early 80s 
[3, ■ These early purely gluonic calculations have been 
improved steadily by including contributions from dy- 
namical quark degrees of freedom (for recent reviews see 
In the most recent calculations the equation of 
state has been evaluated for 2+1 flavor QCD, i. e. in 
QCD with one strange quark and two light (u, d) quarks 
using various improved staggered fermion actions 
The most extensive calculations of the EoS have been 
performed with p4 and asqtad staggered fermion formu- 
lations on lattices with temporal extent N T = 4, 6 0, [1| 
and 8 @. These actions improve both the flavor sym- 
metry of the staggered fermions as well as the quark 
dispersion relations. The latter insures that thermody- 
namic observables are 0(a 2 ) improved at high tempera- 
tures and thus have only a small cut-off dependence in 
this regime. The stout-link action, which has been used 
for the calculation of the EoS on lattices with temporal 
extent N T =4, 6 [6j, only improves the flavor symme- 
try of the staggered fermions and therefore has the same 
large discretization errors at high temperatures as the 
standard staggered fermion formulation. Indeed, the cal- 
culations of the EoS with the stout-link action on lattices 
with temporal extent N T = 4 and 6 reflect the expected 
strong cutoff dependence above the transition tempera- 
ture |||. 

While at high temperatures the masses of the rele- 
vant degrees of freedom, quarks and gluons, are small 
compared to the temperature scale, this is not the case 
at low temperatures and in the transition region. One 



thus may expect that at these temperatures thermody- 
namic observables are more sensitive to the quark masses, 
which control the mass of the light pscudo-scalars and 
eventually are responsible for the occurrence of a true 
phase transition in the chiral limit. Calculations with 
p4 and asqtad actions have so-far been performed us- 
ing light quark masses {m{) which are one tenth of the 
strange quark mass (m s ) and correspond to a pseudo- 
scalar Goldstone mass 1 of 220 and 260 MeV respectively 
Q . The calculations with the stout-link action have been 
performed at the physical value of the light quark mass. 

The purpose of this paper is to investigate the quark 
mass dependence of the EoS by calculating it with the p4- 
action for physical values of the (degenerate) light quark 
masses. The calculational procedure used in this work 
closely follows that used in our previous calculations at 
rhi = 0.1rh s 0]. The paper is organized as follows. In 
the next section we discuss the parameters and some 
technical details of the numerical calculations, including 
the choice of the quark masses and the determination 
of the lattice spacing that fixes the temperature scale, 
aT = 1/N T . In section III we show the QCD equation 
of state and compare it with the resonance gas model. 
In section IV we discuss chiral and deconfining aspects 



1 In calculations with staggered fermions flavor symmetry is bro- 
ken at non-vanishing values of the lattice spacing a. As a con- 
sequence only one of the pseudo-scalar mesons has a light mass 
that is proportional to y'mT and vanishes in the chiral limit at 
fixed a > 0. Full chiral symmetry with the correct Goldstone 
pion multiplet is recovered only for a — > 0. For an estimate of 
the remaining flavor symmetry violations in spectrum calcula- 
tions with the p4-action see lldll . 
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of the QCD transition. Finally section V contains our 
conclusions. In the appendix we give further numerical 
details of our calculations. 



II. ZERO TEMPERATURE CALCULATIONS 

We have performed calculations with the p4-action for 
several values of the gauge coupling f3 = 6/g 2 in the re- 
gion of the finite temperature crossover for lattices with 
temporal extent N T = 8. The values of the quark masses 
and gauge coupling are given in Table |TJ Simulations 
were performed with the Rational Hybrid Monte-Carlo 
(RHMC) algorithm with Hasenbush preconditioning and 
different time steps in the molecular dynamic (MD) evo- 
lution for gauge fields, strange quarks and light quarks. 
We used 5 time steps for strange quark updates per each 
light quark update. The gauge fields were updated 10 
times more frequently than the strange quarks. The 
length of the MD trajectory was 0.5 and the acceptance 
rate about 70%. The zero temperature calculations were 
performed on 32 4 lattices. 

We would like to perform calculations of thermo- 
dynamic quantities keeping the physical values of the 
strange quark mass as well as the light quark masses 
fixed. As a starting point for our calculations we used the 
same fine-tuned values of the bare strange quark mass as 
in 0] but bare light quark masses which are twice smaller 
than in Ref. 0, i.e. the ratio of strange to light quark 
mass was chosen to be h = m s /mi = 20. 

At each beta value we then have determined the static 
quark potential. The static potential V q q{r) and the cor- 
responding scales tq and r\ defined as 



,2 dV q q{r) 
dr 



1.65 , n = 
1.0 , n = 1 



(1) 



were calculated from smeared Wilson loops as described 
in Ref. [ll|. Here we used 30 levels of APE smearing 
of the spatial links to get good signals for the Wilson 
loop expectation values. The values of ro/a and ri/a are 
given in Table HI The comparison with previous results 
given in Table I of Q reveals that the new values of the 
scale parameters in lattice units differ from the old ones 
obtained at h = 10 by less than 1 %, with deviations 
which are not systematic. Thus, the smaller values of the 
light quark masses do not change the lattice spacing in 
physical units beyond the errors, which are about 0.5% 
(c.f. Table QJ. As the consequence the values of the 
temperature in units of ro have errors of about 0.5%. 
To give the temperature in physical units, e.g. MeV, 
we use the value ro = 0.469(7)fm. Thus the temperature 
scale in physical units has an overall uncertainty of about 
1.5%. However, since the same uncertainty is present in 
the calculations performed at twice larger quark mass it 
does not show up in the comparison of thermodynamic 
quantities obtained at two different quark masses. 

To remove the additive divergent constant cq in the po- 
tential, following Ref. Q we normalized it to the string 



form of the static quark potential, V^tring (?*) = — 7r/12r + 
ar, at distance r = 1.5ro. This amounts to rcnormalizc 
the temporal gauge links with z(f3) = exp(co/2). The re- 
sulting multiplicative renormalization factors, z(/3), are 
also given in Table U and will be used to define renormal- 
ized Polyakov loops that are discussed in section IV. 

Having determined the scale we extracted pseudo- 
scalar meson masses using wall sources in the calculation 
of meson propagators. It turned out that the rj S s mass in 
lattice units is, with 1 -2 % accuracy, the same as in 
Thus, a re-adjustmcnt of the line of constant physical rj sS 
mass corresponding to our new and smaller light quark 
masses was not necessary. In fact, in the present calcula- 
tions we find that our quark mass values define a line of 
constant physics characterized by the following relations 



ro • m T = 0.371(3) 

ro ■ m VsS 



ro • mjf = 1.158(5), 
1.578(7) . (2) 



Using ro = 0.469 fm, as determined in Ref.[l2j, we get 
m w = 154 MeV, m K = 486 MeV and 2 = 663 MeV. 
This means that both the light quark masses and the 
strange quark mass are very close to their physical values. 
Furthermore, in the entire parameter range covered by 
our thermodynamic calculations deviations of the meson 
masses from the above values are less than 3%. 



III. CALCULATION OF THE 
THERMODYNAMIC QUANTITIES 

The calculation of the EoS starts with the evalua- 
tion of the trace anomaly, i.e. the trace of the energy- 
momentum tensor 8 W (T). It is related to the temper- 
ature derivative of the pressure through thermodynamic 
identities, 



6 MM (T) _ e - 3p _ T d (p_^ 



2^4 



T 4 



dT \T 4 J 



(3) 



The trace anomaly can be expressed in terms of the ex- 
pectation values of quark condensates and the gluon ac- 
tion density 



e MM (r) _ gro e^(T) 

Rp {( SG )o-(s G ) T }N? 



9£ M (r) 



J--L 



(4) 
(5) 



e£ M (r) 



-RfiRm [2m, «^)j, - W>Z,r) 
+m s {(W>).,o - (V^kr)] V 4 . (6) 



2 A physical value for the r] a s mass can be obtained from the 
relation m„. s = ^2m 2 K — ra\ = 686 MeV. 
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rh s 


#traj. 


ra/a 


ri/a 




rriKCL 


m VsS a 




3.430 


0.0370 


2610 


2.6698(101) 


1.8378(184) 


0.1415(2) 


0.4423(3) 


0.5992(4) 


1.5124(90) 


3.460 


0.0313 


3220 


2.9377(90) 


2.0000(202) 


0.1274(4) 


0.3994(4) 


0.5427(2) 


1.5219(71) 


3.490 


0.0290 


3460 


3.2188(70) 


2.2115(101) 


0.1213(15) 


0.3775(5) 


0.5117(10) 


1.5277(26) 


3.500 


0.0253 


3030 


3.3251(114) 


2.2550(159) 


0.1110(5) 


0.3491(4) 


0.4743(3) 


1.5292(50) 


3.510 


0.0260 


3040 


3.4021(141) 


2.3116(105) 


0.1121(10) 


0.3510(9) 


0.4770(6) 


1.5289(55) 


3.520 


0.0240 


2980 


3.5167(68) 


2.4050(106) 


0.1056(5) 


0.3322(6) 


0.4525(3) 


1.5303(20) 


3.530 


0.0240 


2450 


3.5927(101) 


2.4723(81) 


0.1045(14) 


0.3300(17) 


0.4500(11) 


1.5291(48) 


3.540 


0.0240 


3360 


3.6802(138) 


2.5090(94) 


0.1040(4) 


0.3251(13) 


0.4432(11) 


1.5272(41) 


3.545 


0.0215 


3090 


3.7513(67) 


2.5840(75) 


0.0947(7) 


0.3061(3) 


0.4178(2) 


1.5298(23) 


3.560 


0.0205 


3010 


3.8790(67) 


2.6732(81) 


0.0935(12) 


0.2941(13) 


0.4024(11) 


1.5275(18) 


3.585 


0.0192 


1400 


4.1501(118) 


2.8503(106) 


0.1049(6) 


0.3258(17) 


0.4434(11) 


1.5275(66) 


3.600 


0.0192 


4080 


4.3033(181) 


2.9351(73) 


0.0888(20) 


0.2731(6) 


0.3722(7) 


1.5277(41) 


3.630 


0.0170 


2000 


4.6853(349) 


3.1858(167) 


0.0773(24) 


0.2442(5) 


0.3343(9) 


1.5293(82) 


3.660 


0.0170 


2850 


4.8497(232) 


3.3368(138) 


0.0757(7) 


0.2364(12) 


0.3246(8) 


1.5175(50) 



TABLE I: The parameters of the zero temperature calculations, the values of tq and r\ and the pseudo-scalar meson masses. 
The last column shows the normalization constant for the static potential. 



where the zero temperature expectation values are sub- 
tracted to render the trace anomaly UV finite. The ex- 
pectation values of the light and strange quark conden- 
sates and the action are defined as 



1 1 



AN*N X 
q = I, s , x = 0, t 



(TrD- 1 (m q )) a 
0, r , 
(Sg) x , 



(7) 
(8) 



with D{rhq) being the staggered fermion matrix. Fur- 
thermore, Rp and R m are the non-perturbative beta 
function and the mass anomalous dimension 



W) 



R m {P) 



d(3 

da 

dm, (/3) 



m0) d/3 



(9) 
(10) 



Following Ref. [8[ the /3-dependence of the Sommer scale 
is fitted to a renormalization group inspired Ansatz [l3j 



ro _ 
a 

where 



l + e r a 2 (/3) + / r a 4 (/3) 



OrRnffi) (1 + b r a 2 (l3) + c r a 4 {(3) + d r a 6 (/3)) 



exp 



P 
12&o 



66o 

P 



-b 1 /(2b 2 ) 



(11) 



(12) 



is the two-loop beta function for 3-flavor QCD and 
a((3) = i?2(/3)/^2(3.4). From this the non-perturbative 
beta function Rp can be calculated as 



(13) 




140 160 180 200 220 240 260 



FIG. 1: The trace anomaly (e — 3p)/T 4 calculated for the 
physical quark mass and compared with previous calculations 
at larger light quark masses mi = 0.1m s as well as with the 
HRG model which includes all the resonances up to 2.5GeV. 
Also shown are the interpolations of the lattice data. 



Likewise, R m was obtained from a paramctrization of the 
bare quark mass, 



/19/. \ 4 / 9 
rgi^ I iZb o 



V P 



P(P) 



(14) 



with a sixth order rational function P((3) as in Q to 
account for deviations from the leading order scaling re- 
lation. 

We performed finite temperature calculations on 32 3 x 
8 lattices for the 14 parameter sets shown in Tabic fl] The 
number of MD trajectories for each finite temperature 
run is given in Table ITT] Using the expectation values 
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of the quark condensates and the gluonic action density 
as well as the non-perturbative beta functions described 
above we have calculated the trace anomaly. The nu- 
merical results are shown in Figure [1] and are compared 
to the previous calculation at twice larger quark mass 
mi = 0.1m s on N T = 6 lattices Q and N T = 8 lattices 
[?|. The differences between N T — 6 and N T = 8 calcu- 
lations are due to cutoff effects and have been discussed 
inRcf. [|. 

As one can see from the figure the main differences to 
the iV T = 8 results at m; = 0.1m s arise for temperatures 
T<200 MeV. This difference can be understood as result- 
ing from a shift of the transition temperature when the 
light quark mass is lowered to approximately its physical 
value. Based on our previous study on coarser lattices 
we expect that the reduction of the quark mass from 
0.1m s to 0.05m s effectively leads to a shift of several ob- 
servables by a few MeV towards smaller values of the 
temperature [III ■ As discussed further in section IV this 
shift indeed amounts to about 5 MeV. At lower tempera- 
tures it also is expected that the trace anomaly increases 
with decreasing quark masses as hadrons become lighter 
when the quark mass is decreased. While a tendency for 
such an increase may be indicated by the data at the 
lowest two temperatures reached in our calculation, this 
effect is certainly not significant within the current sta- 
tistical accuracy. The observed consistency within errors 
between the trace anomaly calculated at mi = 0.05m s 
and at mi = 0.1 m s might possibly also hint towards a 
distortion of the hadron spectrum due to finite lattice 
spacing effects. In fact, it is known that for improved 
staggered fermions the ground state hadron masses ap- 
proach their continuum limit from above [Til . Irs] ]. For 
temperatures T > 200 MeV there is no visible depen- 
dence of (e — 3p)/T on the quark mass. This is presum- 
ably due to the fact that hadronic degrees of freedom are 
no longer relevant in this temperature domain; the rele- 
vant degrees of freedom have thermal masses of the order 
of the temperature and are insensitive to the light quark 
masses already for mi = 0.1m s . 

At temperatures below the transition temperature it 
is expected that thermodynamic quantities are well de- 
scribed by a hadron resonance gas (HRG) model. In 
fact, the freeze-out of hadrons in heavy ion experiments 
takes place in the transition region and the observed par- 
ticle abundances are well described by the HRG model 
[HI, [TtJ ■ Therefore in Figure Q] we also show the predic- 
tion of the HRG model, which includes all the known 
resonances up to the mass M max = 2.5 GeV. The lattice 
data for e — Zp are below the HRG prediction although 
the deviations from it are smaller compared to the results 
obtained at m; = 0.1m s . We mention again the present 
statistical accuracy and the possibility of discretization 
effects in the hadron spectrum. In particular, due to taste 
breaking of staggered fermions pseudo-scalar mesons are 
not degenerate at finite lattice spacing, therefore their 
contribution to thermodynamic quantities maybe sup- 
pressed. 



From the trace anomaly the pressure and thus other 
thermodynamic quantities can be calculated by perform- 
ing the integration over the temperature 



— --Tf = ] Ta dT ¥-^ {T) - (15) 



Here To is an arbitrary temperature value that is usu- 
ally chosen in the low temperature regime where the 
pressure and other thermodynamical quantities are sup- 
pressed exponentially by Boltzmann factors associated 
with the lightest hadronic states, i.e. the pions. Energy 
e and entropy (sT = (p + e)) densities are then obtained 
by combining results for p/T 4 and (e — 3p)/T i . 

To perform the integration numerically a reliable in- 
terpolation of the lattice data on the trace anomaly is 
needed. These interpolations are shown in Figure Q] as 
curves. In the region T < 175 MeV, the curves corre- 
spond to exponential fits. Therefore our parametrization 
ensures that the pressure is zero at T = To = 0. How- 
ever, the fits give very small values of the trace anomaly 
already at temperatures around lOOMeV, therefore we 
could also use To = lOOMeV as lower integration limit. 
Above T = 175 MeV, we divide the data into several 
intervals and perform quadratic interpolations. In each 
interval, these quadratic fits have been adjusted to match 
the value and slope at the boundary with the previous 
interval. These interpolating curves are then used to cal- 
culate the pressure and other thermodynamic quantities 
using Eqs. §3§ and (|T5|) . The numerical results for the 
pressure and energy density are shown in Figure At 
temperatures between 170 MeV and 200 MeV, energy 
density and pressure are slightly larger than in earlier 
calculations performed at mi = 0.1to s For temper- 
atures below T < 170 MeV the pressure is almost the 
same. 

Because the discretization errors increase as the tem- 
perature decreases it is interesting to consider other 
choices for the normalization point To. At sufficiently low 
temperatures the HRG model provides a fair estimate for 
the pressure. In particular, at T = 100 MeV the pressure 
is not very sensitive on how many resonances above 1.5 
GeV are included in the HRG model. Therefore we also 
calculated the pressure taking the lower integration limit 
to be T = 100 MeV with p(T )/T 4 = 0.265 calculated in 
the HRG model. The difference in the pressure and en- 
ergy density calculated this way and using the standard 
procedure, where p(To) = 0, can be used as an estimate 
of the systematic error. It is shown in Figure [5] as the 
horizontal band in upper right corner. The thickness of 
the band indicates the expected size of the systematic 
error. The systematic error estimated this way is about 
8% in the energy density at the highest temperature of 
T ~ 260MeV, and about 13% for the pressure. 
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In the opposite limit of zero quark mass one expects a 
chiral transition and the corresponding order parameter 
is the quark condensate defined in section III. For a gen- 
uine phase transition, i.e. in the chiral limit the quark 
condensate vanishes at the critical temperature T c . How- 
ever, we expect that even for the crossover at finite quark 
mass the light quark condensate rapidly drops in the 
transition region, indicating an approximate restoration 
of the chiral symmetry. At non-vanishing quark mass 
the quark condensate needs additive and multiplicative 
renormalization. Therefore, following Ref. [|| we intro- 
duce the so-called subtracted chiral condensate 



140 160 180 200 220 240 260 



(#>i, T -^W)v 



(18) 



FIG. 2: Energy density and three times the pressure at the 
physical value of the light quark mass and compared with pre- 
vious calculations performed at mi — 0.1m s . The horizontal 
band shows the expected uncertainty in the energy density 
due to the choice of the lower integration limit (see text). 



IV. DECONFINEMENT AND CHIRAL 
ASPECTS OF THE QCD TRANSITION 

In this section we are going to discuss the deconfine- 
ment and chiral aspects of the QCD transition in terms 
of renormalized Polyakov loop, subtracted chiral conden- 
sate and strangeness susceptibility. The QCD transition 
in terms of these quantities has been studied previously 
inRefs. |iliJ[ll|. 

In the previous section we have seen that the energy 
density shows a rapid rise in the temperature interval 
T = (170 - 200)McV. This is usually interpreted to be 
due to deconfincment, i.e. liberation of many new de- 
grees of freedom. For sufficiently large quark mass this 
transition is known to be a first order transition (see e.g. 
Ref. [HI). In the limit of infinitely large quark mass the 
order parameter for the deconfinement phase transition 
is the Polyakov loop. After renormalization it can be re- 
lated to the free energy of a static quark anti-quark pair 
F OQ (T) at infinite separation [20H23 



L ren (T) =cxp(-F 0O (r)/(2T)). 



(16) 



A rapid change in this quantity is indicative for deconfinc- 
ment also in the presence of light quarks. The renormal- 
ized Polyakov loop is calculated from the bare Polyakov 
loop by multiplying it by the renormalization constant 
z(/3) given in Table HI 



L ren (T) = z((3) N ^ L bare {f3) = 



(17) 



In the above formula Uq(xq, x) denotes the temporal link 
variables. Note that after performing the ensemble aver- 
age Lb are (f3) is independent of the space coordinate x. 



Subtraction of the strange quark condensate multiplied 
by the ratio of the light to strange quark mass removes 
the quadratic divergence proportional to the quark mass. 

In Figure [3] we show the renormalized Polyakov loop 
and the subtracted chiral condensate A; jS and com- 
pare with previous calculations performed at light quark 
masses equal to one tenth of the strange quark mass 0] . 
The renormalized Polyakov loop rises in the tempera- 
ture interval T = (170 — 200) McV where we also see 
the rapid increase of the energy density. At the same 
time the subtracted chiral condensate rapidly drops in 
the transition region, indicating that the approximate 
restoration of the chiral symmetry happens in the same 
temperature interval as deconfinement. Compared to the 
calculation performed at light quark masses equal to one 
tenth of the strange quark mass we see a shift of the 
transition region by roughly 5 MeV. We note that such 
a shift arises differently in different obscrvables. In the 
case of the subtracted chiral condensate, for instance, a 
major ingredient to the 'shift' is the fact, that at fixed 
temperature the condensate in the transition region is 
strongly quark mass dependent and drops proportional 
to yjraijm s [23|. 

The fluctuation of strangeness is also indicative of de- 
confincment. It can be defined as the second derivative of 
the free energy density with respect to the strange quark 
chemical potential 



Xs(T) 



1 d 2 \nZ{T^ s ) 



T 3 V 



dfi 2 s 



(19) 



At low temperatures strangeness is carried by massive 
hadrons and therefore strangeness fluctuations are sup- 
pressed. At high temperatures strangeness is carried by 
quarks and the effect of the strange quark mass is small. 
Therefore strangeness fluctuations are not suppressed at 
high temperatures. As discussed in Ref. Q strangeness 
fluctuations behave like the energy density in the transi- 
tion region, i.e. they rapidly rise in a narrow temperature 
interval. In Fig. [4] we show the strangeness fluctuations 
calculated at m; = 0.05m s and compare them with pre- 
vious calculations performed at m; = 0.1m s Q. In the 
bottom figure we also show the strangeness fluctuation 
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FIG. 3: The renormalized Polyakov loop (top) and the sub- 
tracted chiral condensate (bottom) as function of the temper- 
ature calculated at m; = 0.05m s and at 0.1m s . 
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FIG. 4: Strangeness fluctuations as function of the tempera- 
ture calculated at mi = 0.05m s and at 0.1m 3 . In the bottom 
figure the numerical data for mi = 0.1m s have been shifted 
by 5 MeV. 



for mi = 0.1m s with a 5 MeV shift of the temperature 
scale. As one can see this shift accounts for most of the 
quark mass dependence of the strangeness fluctuations. 
This is consistent with the conclusion obtained from the 
quark mass dependence of other thermodynamic observ- 
ables. 



V. CONCLUSION 

We have calculated the EoS, renormalized Polyakov 
loop, subtracted chiral condensate and strangeness fluc- 
tuations in (2+l)-flavor QCD in the crossover region from 
low to high temperatures using the improved p4 stag- 
gered fermion formulation on lattices with temporal ex- 
tent N T = 8 at physical values of the light and strange 
quark masses. We found that thermodynamic quantities 
below the deconfinement transition are larger compared 
to the previous calculations performed at twice larger 
quark mass but fall below the resonance gas model re- 
sult. The differences in the thermodynamic quantities 



calculated at m; = 0.05to s and m; = 0.1m s can be well 
understood in terms of the shift of the transition tem- 
peratures towards smaller values when the quark mass is 
decreased. This conclusion is also supported by the cal- 
culation of renormalized Polyakov loop, subtracted chiral 
condensate and strangeness fluctuations. No additional 
enhancement of the pressure and the energy density is 
seen at low temperatures. This and the deviation from 
the resonance gas model may be a cutoff effect due to 
taste violations. However, better statistical accuracy and 
calculations at smaller lattice spacing are needed to quan- 
tify this assertion. At temperatures above 200 MeV no 
quark mass dependence is seen in the equation of state. 
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Appendix 

In this appendix we present some numerical details of 
our calculations. In Table|ll]we give the expectation val- 
ues of the gauge action and quark condensates calculated 
on 32 4 (zero temperature) and 32 3 x 8 (finite tempera- 
ture) lattices. In Table UTTl we present the numerical val- 
ues of the trace anomaly, pressure, energy density, bare 
Polyakov loop and strangeness fluctuations. 



[1] J. Engels, F. Karsch, H. Satz and I. Montvay, Phys. Lett. 

B 101, 89 (1981) and Nucl. Phys. B 205, 545 (1982). 
[2] J. B. Kogut, H. Matsuoka, M. Stone, H. W. Wyld, 

S. H. Shenker, J. Shigemitsu and D. K. Sinclair, Phys. 

Rev. Lett. 51, 869 (1983); 

T. Celik, J. Engels and H. Satz, Phys. Lett. B 129, 323 
(1983). 

[3] C. E. DeTar, PoS LAT2008, 001 (2008). 

[4] P. Petreczky, Nucl. Phys. Proc. Suppl. 140, 78 (2005). 

[5] P. Petreczky, Nucl. Phys. A 830, 11C (2009). 

[6] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, JHEP 
0601, 089 (2006). 

[7] C. Bernard et al, Phys. Rev. D 75, 094505 (2007). 

[8] M. Cheng et al, Phys. Rev. D 77, 014511 (2008). 

[9] A. Bazavov et al, Phys. Rev. D 80, 014504 (2009). 
[10] M. Cheng et al, Eur. Phys. J. C 51, 875 (2007). 
[11] M. Cheng et al, Phys. Rev. D 74, 054507 (2006). 
[12] A. Gray et al, Phys. Rev. D 72, 094507 (2005). 
[13] C. Allton, Nucl. Phys. B [Proc. Suppl.] 53, 867 (1997). 
[14] C. W. Bernard et al, Phys. Rev. D 64, 054506 (2001). 



[15] C. Aubin et al, Phys. Rev. D 70, 094505 (2004). 
[16] J. Cleymans and K. Redlich, Phys. Rev. C 60, 054908 
(1999). 

[17] A. Andronic, P. Braun-Munzinger and J. Stachel, Nucl. 

Phys. A 772, 167 (2006). 
[18] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. 

Lett. B 643, 46 (2006); Y. Aoki, S. Borsanyi, S. Durr, 

Z. Fodor, S. D. Katz, S. Krieg and K. K. Szabo, JHEP 

0906, 088 (2009). 
[19] F. Karsch, C. Schmidt and S. Stickan, Comput. Phys. 

Commun. 147, 451 (2002). 
[20] L.D.McLerran and B.Svetitsky, Phys. Rev. D 24, 450 

(1981). 

[21] O. Kaczmarek, F. Karsch, P. Petreczky and F. Zantow, 

Phys. Lett. B 543, 41 (2002) 
[22] P. Petreczky and K. Petrov, Phys. Rev. D 70, 054503 

(2004); O. Kaczmarek and F. Zantow, Phys. Rev. D 71, 

114510 (2005). 
[23] S. Ejiri et al, Phys. Rev. D 80, 094505 (2009). 



8 






rh s 


T [MeV] 


# traj. 


(sg)o 


(sg)t 


('<H>)i,o 


(lH>)l,T 


{lM>)s,0 


(ipi>)s,r 


3.4300 


0.0370 


139 


7950 


4.10978(19) 


4.10922(15) 


0.06996(19) 


0.06788(12) 


0.14253(12) 


0.1419(8) 


3.4600 


0.0313 


154 


15290 


4.04326(16) 


4.04242(17) 


0.05263(14) 


0.04932(9) 


0.11670(8) 


0.1157(8) 


3.4900 


0.0290 


170 


33710 


3.98357(19) 


3.98215(11) 


0.04018(13) 


0.03453(11) 


0.10023(10) 


0.0984(7) 


3.5000 


0.0253 


175 


32520 


3.96333(14) 


3.96161(12) 


0.03538(12) 


0.02783(18) 


0.08873(9) 


0.0862(9) 


3.5100 


0.0260 


180 


30050 


3.94573(19) 


3.94322(7) 


0.03294(13) 


0.02339(20) 


0.08733(10) 


0.0839(7) 


3.5200 


0.0240 


185 


39990 


3.92713(9) 


3.92431(10) 


0.02956(6) 


0.01718(17) 


0.08003(5) 


0.0756(8) 


3.5300 


0.0240 


191 


70230 


3.90994(11) 


3.90651(6) 


0.02715(9) 


0.01250(16) 


0.07735(7) 


0.0721(7) 


3.5400 


0.0240 


196 


56740 


3.89333(10) 


3.88947(6) 


0.02518(13) 


0.00929(8) 


0.07507(8) 


0.0691(5) 


3.5450 


0.0215 


198 


20700 


3.88427(6) 


3.88048(9) 


0.02328(6) 


0.00681(12) 


0.06848(3) 


0.0618(9) 


3.5600 


0.0205 


206 


11660 


3.85944(14) 


3.85579(7) 


0.02039(11) 


0.00459(6) 


0.06315(8) 


0.0559(8) 


3.5850 


0.0192 


219 


17250 


3.81987(14) 


3.81650(10) 


0.01638(8) 


0.00297(2) 


0.05583(7) 


0.0479(8) 


3.6000 


0.0192 


227 


4260 


3.79720(7) 


3.79408(17) 


0.01478(7) 


0.00272(2) 


0.05359(4) 


0.0460(9) 


3.6300 


0.0170 


243 


10810 


3.75307(9) 


3.75040(9) 


0.01141(8) 


0.00202(7) 


0.04524(6) 


0.0378(4) 


3.6600 


0.0170 


259 


45050 


3.71067(8) 


3.70856(3) 


0.00916(11) 


0.00184(2) 


0.04202(5) 


0.0355(2) 



TABLE II: The expectation values of the gauge action, light and strange quark condensates at zero and finite temperatures. 
Also shown are the temperature values for each value of gauge coupling p obtained from ro = 0.469fm. We also give the number 
of MD trajectories for each run 



T [MeV] 


(e - 3p)/T 4 


P/T 4 


e/T 4 




Xs/T 2 


139 


0.89(30) 


0.067 


1.09(30) 


0.00087(6) 




154 


1.34(30) 


0.161 


1.82(30) 


0.0015(5) 


0.062(4) 


170 


2.34(30) 


0.317 


3.29(30) 


0.0028(6) 


0.150(9) 


175 


2.87(27) 


0.393 


4.05(27) 


0.0036(9) 


0.186(10) 


180 


4.23(29) 


0.495 


5.71(29) 


0.0046(9) 


0.261(10) 


185 


4.87(19) 


0.621 


6.73(39) 


0.0057(6) 


0.335(6) 


191 


5.99(19) 


0.768 


8.30(19) 


0.0069(7) 


0.392(10) 


196 


6.84(18) 


0.929 


9.63(18) 


0.0082(4) 


0.462(5) 


198 


6.78(18) 


1.013 


9.81(18) 


0.0089(9) 


0.492(7) 


206 


6.73(27) 


1.274 


10.55(27) 


0.0104(11) 


0.599(11) 


219 


6.48(31) 


1.703 


11.59(31) 


0.0135(17) 


0.750(10) 


227 


6.15(34) 


1.926 


11.93(34) 


0.0146(24) 


0.798(8) 


243 


5.49(26) 


2.296 


12.38(26) 


0.0172(11) 


0.841(7) 


259 


4.58(18) 


2.602 


12.38(18) 


0.0199(5) 


0.874(3) 



TABLE III: The numerical values of the trace anomaly, pressure, energy density, bare Polyakov loop and strangeness fluctuations 
at each temperature. 



